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Abstract 

The spreading (propagation) of diseases, viruses, and disasters such as power blackout through a 
huge-scale and complex network is one of the most concerned issues today. In this paper, we study 
the control of such spreading in a nonlinear spreading model of small-world networks. We found 
that the short-cut adding probability p in the N-W model jTj of small-world networks determines the 
Hopf bifurcation and other bifurcating behaviors in the proposed model. We further show a control 
technique that stabilize a periodic spreading behavior onto a stable equilibrium over the proposed 
model of small- world networks. 
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1 Introduction 

Recently, Watts and Strogatz proposed a small- world model ^I], named W-S model, which consists of a 
rewired regular lattice with a very small fraction of long-range connections. Newman and Watts (T) then 
made a slight modification on the rewiring process of the W-S model, named N-W model, by means of 
adding linkages between pairs of randomly chosen nodes with a very small probability < p ^ 1. When 
selecting a very small probability < p ^ 1 in these models, the obtained small-world networks have 
large clustering and small average distances, which match the familiar small-world feature discovered in 
many real-life large-scale networks |31 |H1 El EI ■ 

Since the discovery of the small-world phenomenon, many interesting results on the analysis, control 
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and applications of small- world networks have been published El EH El- In particular, Moukarzel 
|S] studied a linear disease spreading model of small-world networks, which was further extended to a 
nonlinear model with frictions and time-delays where chaos and bifurcation can emerge (21 . In this 
paper, we concern with the Hopf bifurcation phenomenon during the spreading process in small-world 
networks, and its control problem, by utilizing the short-cut adding probability p as the means to stabilize 
the bifurcating behaviors. 

The remainder of this paper is organized as follows: Section 2 describes the nonlinear spreading model, 
where the existence of Hopf bifurcation in this model is discussed. Section 3 proposes a method of of 
controlling the bifurcating behaviors using the short-cut adding probability p. The effectiveness of this 
control technique is visualized via a numerical example, also in this section. Finally, Section 4 concludes 
the paper. 



2 Nonlinear spreading model and Hopf bifurcation 

Assume that a disease (virus, power blackout) spreads with a constant radial velocity v = I from an 
original infection site of a network. The infected volume V{t) grows according to the following nonlinear 
differential equation: 



dt 



1 + 2pV{t -S)~ + 2p)V^{t - 5) 



(1) 



where 5 is the time-delay during the spreading. The nonlinear friction term consists of two parts: the 
former comes from the nonlinear interaction within the regular lattice, while the latter is the nonlinear 
interaction coming from the newly added links, which is dependent on the probability p in the N-W 
small- world model |21. And /i > is a measure of such nonlinear frictions. 

It should be noted that Eq. |Q contains as special cases all the existing linear models of Newman and 
Watts [71 and Moukarzel 0. When /i = and (5 = 0, obviously the model produces exponential growth 
without limit. Besides, the nonlinear model of Yang |T21 hides the effect of the probability p on bifurcation 
behaviors in spreading, while our model is more general and uncovers this effect. Therefore, model 
(jH) is more suitable for studying the spreading phenomenon in small- world networks. 

Denote s{t) = V{t) — V* , where V* is the equilibrium of Eq. 1^1 with V{t) > and is given by 



V* 



P+Vp^+ + 2p) 
nil + 2p) 



(2) 



To study the stability of the spreading behaviors in model , we select as the parameter of bifurcation 
for analysis. We have the following theoretical results for the existence of Hopf bifurcation in model (Q, 
including the directions, stabilities and periods of the bifurcating periodic solutions. 

Theorem 1 If 5 < j^, then when the positive parameter fi passes through the critical value /xq = 

1652 (l+2p) 



TT 16(5 u n 

— — there is a Hopf bifurcation in model 0) at its equilibrium V* 



Proof: See reference l^. 



To state the next theorem, define 
B = 



g2o = 2^i{l + 2p)B, 



1 - 2^p^ + + 2p)5e-*"«* 
gn = -2^(1 + 2p)B, = 2Ai(l + 2p)B, 



321 = -2/i(l + 2p)Bi 



,511 +.gii + 2^(1 + 2p) -.920 - .902 + 2Ai(l + 2p) 



2Vp2 + ^(l + 2p) 2iLJo + 2v/p2 +/i(l + 2p) 
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Figure 1: The evolution trend of bifurcation parameters Hoip), l^2{p), T2{p), /32(p), with < p < 0.2 in 
model (small- world networks), where S = ^. 
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a'(0) = ^ [Re A] 



(1 + 2p)S 



M2 



Re Ci(0) 
a'(0) ' 



cj'(O) = — [Im A] 



T2 



-Im Ci(0) + M2^'(0) 

0^0 



P2 



_ 16(5(1 2p) 
.0 " ^(4 + 7r2) ' 

2Re Ci(0). 



Theorem 2 Parameter /i2 determines the direction of the Hopf bifurcation: If fi2 > (< 0), the Hopf 
bifurcation is supercritical (suhcritical) and bifurcating periodic solutions exist for /i > /ig (< /io); pa- 
rameter P2 determines the stability of the bifurcating periodic solutions: the solutions are orbitally stable 
(unstable) if (32 < (> 0); and parameter T2 determines the period of the bifurcating periodic solutions: 
the period increases (decreases) if T2 > (< 0). 



Proof: See reference [2]. 



3 Control of bifurcating behaviors in small-world networks 

As can be seen from Theorems I1I2I the critical value of /io for the occurrence of Hopf bifurcation, and 
those values of /i2, r2, (32 for verifying the stability of bifurcating periodic solutions, are all dependent 
on the probability p, which will be denoted by Hoip), M2(p), T2(p), /92(p) below. 

Generally, small- world networks are described by W-S/N-W models [3^2 with a small probability 
< p <C 1. Figure n shows the values of ^J,a{p), M2(p), T2{p), P2{p) when < p < 0.2, where, interestingly, 
we have found that T2 (p) and (32 (p) sue not sensitive to the change of the probability p. It implies that these 
small- world networks have almost the same bifurcating periodic behaviors, for example, r2(0.01) = 0.2162, 
which is almost the same as r2(0.1) = 0.2135. This means that when /j, increases after passing through 
the critical value /zo(p = 0.01,0.1), the periods of these periodic solutions increase in about the same 
amount. 

Recall that bifurcation control deals with a modification of some bifurcation characteristics of a pa- 
rameterized nonlinear system by a judiciously designed control input Because p determines the 
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Figure 2: Phase plot and waveform plot of a small-world network of model with p — 0.1 (dotted 
line), and p = 0.01 (solid line), both 5 = ^ and /i = 1.45 {p changes a.t t = 100). 



bifurcating values of /io(p)i ^2(^)1 which are monotonically decreasing as p increases with a fixed time- 
delay S. Therefore, if we want to apply a controller to the bifurcating behaviors, p is preferred to be 
controlled. We thus employ the following control strategy: 

dV{t) 



dt 



= 1 + 2u[p)V{t -5)- /x(l 2u{p))V\t - 5), 



(3) 



where 

u{p) ^ p — Ap. (4) 

We may select p as the control parameter to first calculate the bifurcation parameter /xq (p) , and then to 
vary p with Ap so as to obtain a different fJ-o{p). This results in different bifurcating spreading behaviors 
in a small-world network. To this end, we can further stabilize periodic spreading solutions onto desired 
stable equilibria, and vice versa. 

Figure |2] shows an example. At the beginning, t G [0, 100], we set p ~ 0.1, 6 = ^, /i = 1.45 in model 
(^, and then compute /io(O.l) = 1.2938 < fi = 1.45. As expected, a Hopf bifurcation occurs and the 
spreading behavior is a periodic solution. To stabilize this periodic orbit onto a stable equilibrium, we 
vary the probability p to be 0.01 after t — 100. Since /xo(O.Ol) = 1.5315 > /i = 1.45, the periodic 
spreading solution asymptotically converges to the stable equilibrium V* — 0.8291, as shown in Fig. [3 



4 Conclusion and discussion 

In this paper, we have introduced a general nonlinear spreading model of the small-world network type, 
which has a flexible nonlinear interaction effect. Based on the study of Hopf bifurcation in this spreading 
model, a simple bifurcation control example has been simulated, showing that the probability p can be 
used as a control parameter to stabilize a periodic spreading behavior onto a stable equilibrium in the 
small- world network. 

More recently, Motter ^ studied a means of controlling cascade failure in a complex network, and pointed 
out that selective removal of network elements can prevent the cascade from propagating through the 
entire network. Such a selective removal control strategy is actually a complementary approach of the 
technique discussed in this paper, for controlling the spreading in small-world networks. 
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